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Abstract 

We use a Boltzmann equation to determine the magnetoconductivity of quan- 
tum wires. The presence of a confining potential in addition to the magnetic 
field removes the degeneracy of the Landau levels and allows one to associate 
a group velocity with each single-particle state. The distribution function 
describing the occupation of these single-particle states satisfies a Boltzmann 
equation, which may be solved exactly in the case of impurity scattering. In 
the case where the electrons scatter against both phonons and impurities we 
solve numerically - and in certain limits analytically - the integral equation 
for the distribution function, and determine the conductivity as a function of 
temperature and magnetic field. The magnetoconductivity exhibits a maxi- 
mum at a temperature, which depends on the relative strength of the impurity 
and electron-phonon scattering, and shows oscillations when the Fermi energy 
or the magnetic field is varied. 

PACS numbers: 72.10.Di, 72.20.My, 73.20.Dx, 73.61.Ey 
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I. INTRODUCTION 



The possibility of making quasi-one-dimensional quantum wires, due to advances in mi- 
crofabrication technology, has greatly stimulated the interest in the transport properties of 
such low- dimensional systems. A number of recent papers have treated the magnetoconduc- 
tivity of quantum wireall using a variety of different theoretical approaches. The suppression 
of scattering between edge-states, due to the presence of the magnetic field, was discussed 
by several authors@a. The effects of impurity and boundary scattering on the Hall effect in 
quantum wires were treated by Akera and Andoi, starting from a Boltzmann equation. The 
Kubo methodi has been used for considering the effect of the coupling to optical phononsfll 
and the influence of impurity scattering! on the magnetotransttort through quantum wires. 
Momentum- and energy-balance equations have been derived0'0 for a quasi-one-dimensional 
gas of electrons in a magnetic field. Other authors0 have used the Keldysh method to dis- 
cuss magnetotransport in the presence of impurity scattering, while the formalism due to 
Landaueilll and ButtikeiH was employed in Ref. [15] to treat the influence of disorder on the 
Hall effect in quantum wires. 

Transport in strong magnetic fields has traditionally been treated within the Kubo- 
formalism, since this allows one to take fully into account the quantization of the motion of 
an electron in a magnetic field. The resulting diagrammatic expansion of the conductivity 
requires the consideration of both self-energy effects and vertex corrections. It was demon- 
strated in Ref. [5] that it is essential to include vertex corrections in order to determine the 
effect of impurity scattering on the magnetoconductivity of a quantum wire. The physical 
reason for this is that only back-scattering from the impurities contributes to the resistivity. 

In the present paper we use a Boltzmann equation to determine the transport properties 
of quantum wires formed by additional confinement of the two-dimensional electron gas of 
e.g. a GaAs-GaAlAs heteroj unction. The presence of a confining potential in addition to 
the magnetic field removes the degeneracy of the Landau levels and allows one to associate 
a group velocity with each single-particle state. The distribution function describing the 
occupation of these single-particle states satisfies a Boltzmann equation, which we solve 
exactly in the case of impurity scattering. The resulting magnetoconductivity agrees with 
that obtained in Ref. |9| within the Kubo-approach. In the case of electron-phonon scattering 
we solve numerically - and in certain limits analytically - the integral equation for the 
distribution function. The conductivity is determined as a function of temperature and 
magnetic field in the case when both impurities and phonons contribute to the scattering of 
the electrons. Our paper extends previous work discussed above by giving a more complete 
account of the simultaneous influence of electron-impurity and electron-phonon scattering. 
By starting from the Boltzmann equation we include vertex corrections from the outset, 
thus simplifying the formal development in comparison to that involved in the use of the 
Kubo formalism. 

Apart from its simplicity, the use of the Boltzmann equation has the additional advantage 
that it allows one to go beyond linear response theory in a straightforward manner. Although 
the present paper treats in detail only the regime of linear response, our approach may readily 
be generalized to take nonlinear effects into account. 

Since conduction in a quantum wire involves current transport in one dimension, it is 
necessary to consider the role of localization. This was done in Ref. [16|, where it was shown 
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that weak localization in quantum wires is destroyed in a magnetic field greater than a 
critical value B c , where the critical magnetic field B c apart from numerical factors is given 
by B c ~ M/eL^W 2 . Here £ is the elastic scattering length, W the width of the quantum 
wire, while is the phase coherence length. For the parameters considered in Ref. [16] this 
yields a B c somewhat less than 0.1 T. In what follows we shall assume that the magnetic field 
is sufficiently strong that localization effects are always negligible. The transport properties 
of the quantum wires may then be obtained from a Boltzmann equation, in which the 
quantizing effect of the magnetic field is incorporated in the manner described below. We 
shall also disregard Coulomb-blockade effects which have been shown to be of importance 
in some systems0. 

The paper is organized as follows. In section 2 we introduce the single-particle states 
derived from the presence of a confining potential in addition to a homogeneous magnetic 
field. The Boltzmann equation and its solution is discussed in section 3. In section 4 we treat 
scattering due to impurities, while scattering from impurities and acoustical phonons forms 
the subject of section 5, where we consider both deformation-potential and piezoelectric 
coupling. Finally section 6 deals with the case where scattering from impurities and optical 
phonons is important. 



II. THE SINGLE-PARTICLE ENERGY SPECTRUM 

We consider an electron moving in the xy-plane under the influence of a constant mag- 
netic field in the ^-direction. In addition a parabolic confinement potential limits its motion 
in the ^/-direction whereby a wire is formed in the x-direction. The Hamiltonian for such an 
electron is 

ir = -L(p + e A)» + (i) 

m* being the effective (band) mass of the electron. For simplicity we have neglected the 
Zeeman splitting, since we shall be mainly concerned with quantum wires in GaAs-based 
structures, where the Zeeman energy is a few percent of the cyclotron energy, tkv c , to be 
defined below. The electron spin therefore only enters as a factor of 2 in the expression 
for the current density. In the x-direction we impose periodic boundary conditions, and for 
the vector-potential A we use the Landau-gauge, A = B(—y, 0, 0). Consequently the eigen- 
functions of the Hamiltonian are products of plane waves and harmonic-oscillator functions 
involving Hermite polynomials H n , 

Vnk(x,y) = -^=e ikx 1 



with L being the length of the quantum wire. The center coordinate, y^, and the charac- 
teristic length, £h, of the harmonic oscillator is specified below. The wave function in the 
^-direction is taken to be a delta function. 
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In the case where K is zero and only the magnetic field is present the natural system of 
units is based on the cyclotron frequency uo c and the magnetic length £ c , given by 

u c = —, e c = -s- 3 

m* eB 

When K is different from zero, we define the quantity 7, which describes the relative strength 
of the confinement potential, by 

7 = (4) 

and the natural units are now based on the hybrid frequency, u>h, and the hybrid character- 
istic length, £h, given by 

u h = {l + ^u) c , 4 = (1 + 7)^4. (5) 
By inserting the function (0) in the Schrodinger-equation 

H^ nk = E nk y nk , (6) 

one finds that @ is satisfied, provided 

e nk = hu h (n + -J + n = 0,1,2,..., (7) 

where m is a renormalized mass defined by 

1 + 7 * , a s 
m = m . (8) 

7 

The center of the harmonic-oscillator wavefunctions is given by 

^ = -^-1 = (l + 7)~^4 (9) 

The confinement potential Ky 2 /2 removes the degeneracy of the Landau-levels through 
the fc-dependence exhibited in (|7|). Note that the dispersion relation is that of a free particle 
with a renormalized mass m. When the strength of the confinement potential goes to zero, 
the renormalized mass m becomes infinite (corresponding to degenerate Landau-levels), 
while in the opposite limit (7 — > 00), the renormalized mass becomes equal to the effective 
mass m*. The fc-dependent energy levels are plotted in Fig. 1, where we have also indicated 
the position of the Fermi level corresponding to a definite electron density. 

In the following we shall describe conduction through the quantum wire on the basis 
of the Boltzmann equation for the distribution function of the excitations with energy e n k- 
The group velocity v of the excitations is seen to be 

^ ^ _ 1 de nk _ hk 
n h dk m 

If the confinement potential is not simply parabolic, the relation between v n k and k be- 
comes more complicated than Eq. (|i~0D, but the method described in the following is still 
applicable, provided one modifies the group velocity accordingly, and incorporates the new 
wavefunctions into the matrix elements appearing in the collision probability. 
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III. THE BOLTZMANN EQUATION 



We consider the distribution function for the excitations specified by the energy e nk , 
with group velocity v nk = hk/m. Unlike the usual semi-classical description of transport in a 
magnetic field, the magnetic field has here been taken into account from the very beginning, 
in defining the excitations with dispersion relation ([?]). The effect of an electric field E in 
the x-direction is, however, included in the usual manner through the acceleration equation 

hk = -eE. (11) 



The Boltzmann equationc3 is 



dfnk . ;df n k ( dfnk\ /1oX 

-w +k sr = \-w)^ < 12 > 

where the right hand side of the equation contains the collision term. Since we are interested 
in the linear response to a static electric field, the time-derivative df / dt is equal to zero, while 
f n k in the second term on the left hand side may be replaced by the equilibrium function /° fe . 
Using Eq. ([TT]) the Boltzmann equation then becomes an inhomogeneous integral equation 
of the form 

( E v nk f° nk (i - a = f%l , as) 



k n T 



with the integral term on the right hand side - the collision term - to be specified below. 

Before we consider specific scattering mechanisms, let us write down the conductivity 
a in terms of the distribution function. In the present work we restrict ourselves to the 
calculation of the diagonal conductivity element a xx = a. The current density j is 

/DO dfc 
— V nkfnk- (14) 
-oo Z7T 

We shall introduce the deviation function tp n k by the definition 

fnk = fnk + ~ fnk^nk- (15) 

Then the current density becomes 

1,00 dk 



J 



/OO rllr 
^ nk fl h {l-f nk )^ nk . (16) 



Since we shall be dealing with distribution functions that change sign upon k — > —k, we may 
restrict k to positive values and work instead with the distribution function as a function 
of energy: ip n (e). In this notation k is then a function of energy, and we denote the k- value 
which solves e = e nk by k n {e). Furthermore, by introducing the function n defined by 

K = -^r^n, (17) 
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we may express the conductivity as 



since dk = de/%v. Note that <fi has dimension of length. The chemical potential entering 
the equilibrium distribution function f° will be assumed to be independent of temperature, 
equal to £p, since the quantum wire under typical experimental conditions is in contact with 
a large reservoir of electrons. 

The Boltzmann equation considered in this paper is an inhomogeneous integral equation 
of the form 



X nk = Hip nk , (19) 

where 

eE 

*n fe = -^W° fe (l-/„° fe ), (20) 



while the integral operator H is defined by 

/ dfnk\ 



9t ; coll 



-Hip nk (21) 



with 

/OO 
—dk'K nn ,(k, k'){ip nk - VVfc')- ( 22 ) 
-oo Z7T 

Here K is an integral kernel to be specified in later sections. 

By introducing the scalar product (A, B) through the definition 

/oo dfc 
—A nk B nk (23) 
-oo Zll 

the conductivity may be written as 

* = %?W)- ( 24 ) 

In cases where we cannot solve explicitly for ip, it is useful to employ a variational principle, 
which yields a lower bound on a by virtue of the Schwarz inequality 

(^H^)(U,HU)>(U,Hij) 2 , (25) 

where U is an arbitrary trial function. Since Hip = X, this gives a lower bound on the 
conductivity 

2k B T (X, Uf 

° - & (u,huy (26) 
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This lower bound will be used to determine the low-temperature conductivity when the 
scattering is due to acoustic phonons, both with and without impurities, and the resulting 
analytic expressions are compared to the results of the numerical solution of the integral 
equation. 

We shall always write our calculated values of the conductivity, as limited by the various 
scattering mechanisms to be considered in subsequent sections, in the form 

2e 2 

*= T l, (27) 

where / has the dimension of a length. As may be expected, in the presence of impurity 
scattering I tends to a finite value as T tends to zero, while it increases exponentially in the 
absence of impurity scattering, when the electrons are scattered only by phonons. 



IV. IMPURITY SCATTERING 

In this section we consider the case of elastic impurity scattering treated recently in 
Ref. H Our aim is to show how the results of Ref. |], which were obtained by use of the 
Kubo-formalism, are derived within the present framework. The collision integral is 

/DO £, 
—dk'w nn ,{k, k')(f nk - f n > k >), (28) 
-oo ZTT 

since the scattering is elastic. According to the Golden Rule we have 

w nn ik,k') = ^-\(nk\V\n'k')\ 2 5{e nk - e n , k ,). (29) 

The square of the matrix-element (nk\V\n'k') for scattering from impurities may be written 

as 

\(nk\V\n'k')\ 2 = J d 2 r J d 2 r'V(r)V(r') 

x (nk\r) (r\n'k f ) (n'k'\r'} (r'\nk). (30) 

We shall perform an ensemble average over the distribution of impurities, corresponding 
to the replacement 

V(r)V(r') -> (V(r)V(r')} = F(r,r'). (31) 

When the impurity potential may be approximated by a delta-function in space, the function 
F is proportional to the delta-function 8{r — r'), corresponding to a "white-noise" model 

F(r,r') = A 2 n imp 5(r-r'). (32) 

Here we have introduced the number of impurities per unit area ni mp together with the 
constant A, which denotes the magnitude of the matrix-element for scattering from a single 
impurity. The dimension of A is that of an energy times an area. 
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In order to distinguish between states of equal energy but opposite sign of k, we introduce 
the "branch-index" s = ±1. 



2m 

sk n (e) = sJ -^-(e - (n + ^)huj h ) 



1 

2' 



(33) 



provided 



(n + -)hoj h < e. 



(34) 



The density of states at the Fermi energy Ep is inversely proportional to the s-independent 
Fermi velocity v n (ep) = ^k n (e). Finally since (p ns is odd in k and hence changes sign with 
s, we may write <p ns = S(j) n . The problem then becomes that of solving the matrix equation 



Bl 



(35) 



and calculating the conductivity according to Eq. (p!8|). The matrix B™/ is positive and 
symmetric with respect to interchange of ns with n's'. Its elements are 



B n's'_}_ A 2 1 1 



em's' 
ns ' 



Here we have defined the following dimensionless quantities 



(36) 



2tt4 / dy\u n (y - y kns )\ \u n >(y-y k )| . 



(37) 



In Fig. 2 we show the result of solving the matrix equation fl35]) and calculating the 
conductivity according to Eq. (|T8| ) as a function of the Fermi energy £ F , for different choices 
of the parameter 7 determining the strength of the confining potential. The pronounced 
structure is due to the opening of new "channels" for scattering each time the Fermi energy 
exceeds (n+l/2)hujh- At finite temperature this structure is somewhat smeared, as evidenced 
by the figure, but still clearly visible. 

When TiiOh/2 < e-p < ?>Tiu} h /2 only the lowest Landau-level is involved, corresponding to 
n = n' = in Eq. (^). Then the solution becomes particularly simple, 4>q = 1/2B®+. By 
means of the Sommerfeld-expansion we obtain the following low-temperature expression for 
the conductivity a = 2e 2 /hB^, 



2e 2 



a 



-lhV2ir 



% 7 



h 2 kl 



A 2 n imp m? 



exp 2^/[l+ 7 ; 



x 



57T 



3r 



1 + 



1+7 ^ f k B T ' 



h 2 l 2 



(38) 



where kp = k 0+ (ep) is the Fermi momentum. The conductivity increases with temperature 
because electrons are excited into states which have a smaller overlap with the states at the 
opposite branch. Consequently the electrons have a smaller back-scattering probability. The 
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result (^) is valid only to lowest order in {k&T ' /TiuJh) 2 - In obtaining our numerical results 
we do not use the Sommerfeld-expansion, but calculate directly the integral over energy 
according to the expression fllB|) . 

As an illustration of the matrix-inversion involved in the solution of the Boltzmann 
equation we shall also give the expression for the zero-temperature conductivity in the case 
when 3huh/2 < £f < 5huh/2. By utilizing the symmetry <ft n+ = —<j) n - the 4x4 matrix- 
equation becomes a 2 x 2 equation which may readily be solved for 0o+ and <f>i+. The 
resulting conductivity o = (2e 2 /h)2((p 0+ + is at T = K given by 

o = ^ 2(2^+ + B\- + + Bt + ) ^ 

h 2(B^Bq+ + Bq+B\+) + + i?o+)(-Bo+ + B\+) 

The result (^) yields the conductivity at zero temperature, when the Fermi energy e-p 
satisfies the condition 3huih/2 < Ep < 5h(Jh/2, in agreement with Ref. |^. The generalization 
to more Landau- levels is straightforward. The B- matrix in Eq. ( p5|) is in general a 2(n + 
1) x 2(n + 1) matrix given by Eq. (|36|), where n is the quantum number associated with 
the highest occupied Landau-level at T = K. 

It is convenient to relate the conductivity calculated in the following two sections to the 
zero-temperature conductivity &2D associated with motion in two dimensions in the absence 
of a magnetic field and a confinement potential. The latter is given by 

0"2D = —Timp = e^DyUimp, (40) 

m* 

where n2D is the electron sheet density and /ii mp is the zero-field, zero-temperature mobility, 

/W = ^ (41) 

m* 

with 

1 m*A 2 n 



imp 



Timp fl 



3 



(42) 



We shall thus express our calculated conductivity in units of (2e 2 /h)k mp , where 

^imp ^imp^c^c (43) 

is a characteristic length which depends on the impurity content as well as the magnetic 
field. 



V. ELECTRON-PHONON SCATTERING 

Next we turn to the consideration of inelastic processes. The present section treats 
the case where the electrons are scattered by acoustical phonons due to the combined 
deformation-potential and piezoelectric coupling, while the following section discusses the 
coupling of electrons to optical phonons. 
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Before specializing to a particular model let us write down the contribution to the collision 
operator from the scattering against phonons. The phonons are assumed to be in thermal 
equilibrium. Then the integral operator H is given by the following expression 

Hl/; nk = -r- \9q\ 2 f°( £ nk)(l ~ f {s n 'k'))h,k'+q x i^nk ~ Ipn'k') 
11 n' k' q 

x{5(e n /^ - e nk + hu q )(l + N°(u q j) + 6(e nlk > - e nk - huj q )N°(uj q )}. (44) 

Here g q is the electron-phonon coupling parameter, which will be specified below, while 
N°(u q ) is the equilibrium phonon distribution. For simplicity we have assumed that the 
phonon frequencies only depend on the magnitude of q. In Eq. (fB|) we have explicitly 
written the factor 8 kyk i +qx discussed in the appendix. Utilizing the relationship 

f(e)(l-f(e±hu> q )) = (f(e)-f(e±?iw q ))(l+N (±u; q )) 1 (45) 

the collision integral (f44|) can be written as 

HI —T\ |2 l/°( g ^fc) - f°{e nk + ahu q )\ 

4sinh 2 (^/2fc B T) 

X (lf) n k -^n'k-q x )K £ n'k-q x ~ £nk ~ (T%U q ), (46) 

with cr assuming the two values +1 and —1 corresponding to phonon absorption and 
emission respectively. 

Our calculations of g q are based on the standard electron-phonon interaction Hamilto- 
nian, which disregarding umklapp processes takes the form 



\ 2Vpu, 



" 52M x (q)p(q)(a q + al q ), (47) 



1 \,q 



where p is the ion mass density, V the normalization volume, p(q) the Fourier-component of 
the electron density operator, aL q a phonon creation operator, and A the polarization index. 
The coupling function M is given by 

Mx(q) = -V(q)q^ x , (48) 

V(q) being the electron-ion potential and £a a unit vector describing the polarization A. For 
a detailed treatment of electron-phonon coupling in semiconductors see e.g. Ref. [19]. The 
values of the GaAs parameters used in our calculations are listed in Table |. 



A. Coupling to acoustical phonons 

In GaAs-heterostructures at low temperatures the electron-phonon scattering is mainly 
due to the combined piezoelectric coupling and the deformation-potential couplingnHc3. 
Below we briefly sketch how these couplings are derived. 

The deformation potential coupling only involves the longitudinal acoustical (LA) pho- 
nons. The coupling to the transverse acoustical (TA) phonons is suppressed by the square 
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of the ratio between the speed of sound and the speed of light. In the long- wavelength limit 
the coupling function M d f for the deformation-potential coupling is written as 



M df =Eq, (49) 

where S, known as the deformation potential, is the zero-wave- vector limit of V{q). 

The electron-ion potential V for the piezoelectric coupling is found from the basic piezo- 
electric equations^. For GaAs (zinc-blende structure) the only non-vanishing independent 
piezoelectric constant is hu (reduced notation), and the coupling function M^ z in this case 
becomes 

M{ z = i2ehu(q x q y £x,z + Wzix,* + q z qx£\,y), (50) 

where = (q/q)i and = m the piezoelectric case A is retained. 

It is noted that M * is real while M pz is imaginary; thus to second order the two terms 
do not interfere, and the absolute square of the total coupling function M ac is given by 

\M ac \ 2 = \M df \ 2 + \M pz \ 2 . (51) 



To obtain a more tractable form of the piezoelectric coupling we perform angular averages 
for the longitudinal and the (two) transverse modes separately and then add the termsu'El. 
While this represents an approximation compared to retaining the full q-dependence of 
the coupling in the collision integral, we expect it to involve only minor quantitative dif- 
ferences. In adding the transverse and longitudinal contribution we must remember the 
different average sound velocities associated with each of the terms, originating in the fac- 
tor 1/ y/UJ^ = l/yfcq in Eq. (|47|) . Expressing the transverse sound velocity as x times the 
longitudinal sound velocity we obtain the following angular average of the absolute square 
of the piezoelectric coupling function M pz : 

\M pz \ 2 =(eh u ) 2 (^ + -^)^P, (52) 



35 x 35 

where we have introduced the constant P, and where it is understood that the only sound 
velocity appearing in the following is the longitudinal one. 

The electron-phonon coupling parameter g q introduced in Eq. ( fH| ) now becomes 

\gUq)\ 2 = - q {^\ 2 + P) \{nk\e^\n'k')\ 2 . (53) 

The matrix element appearing in this expression is treated in detail in the appendix. In 
writing Eq. (|53D we have taken the phonon frequency u q to be given by 

u q = cq. (54) 



11 



B. Coupling to acoustical phonons in the low-temperature limit 



In an analytical study of the low-temperature limit we use the variational principle 
discussed in section 3 above and choose a trial function given by 

U = U nk = sgn(fc). (55) 

Thus U is 1 on the branches corresponding to k being positive, while it is —1 on the branches 
corresponding to k being negative. This choice will lead to an expression for the conductivity 
which in the case of a single occupied Landau-level and in the limit of low temperature agrees 
with a numerical solution of the integral equation. This suggests that the trial function is 
in fact the exact one under these conditions. 

First we evaluate the scalar product (U, X) appearing in the general expression (|26"D for 
the lower bound on the conductivity. By changing the integration variable to the energy e nk 
we get 

eE 

(X,U) = —J2f°(s n0 ), (56) 

where e n o is the k = value of e nk . Because of the symmetry of the integral operator H 
given in Eq. (flop the denominator occurring in Eq. ( ^6|) becomes 

00 dk r°° dk! 



„ J— oo Z7T J -co Z7T , 



xi[l-sgn(A;)sgn(fc')] 2 . (57) 
As discussed in the appendix k' = k — q x . The occurrence of the factor [1 — sgn(fc)sgn(fc — q x )] 



in the integrand of Eq. (]57|) implies that the summation over q x becomes restricted by q x > k 
for k > 0, while q x < k for k < 0. Using the symmetry with respect to reversal of all momenta 
variables we may thus restrict the /c-integral to the interval to oo, which limits q x to the 
region q x > k. If 8 denotes the angle between q and the x-axis, the integration over 6 is 
therefore restricted by 

k 

cos#>-, < q < oo. (58) 

q 

Thus we have to carry out three integrations (over q, k and cos 9) as well as two sums (over 
n and n'). 

Let us consider the simplest case, in which n — n' — 0. If the sound velocity is much less 
than the Fermi-velocity then the energy-conserving delta-functions are 

5(e n > k ,-e nk ±hu} q ) ~ 6( — (q 2 cos 2 9-2kq cos9)). (59) 

2m 

The g-integral is now restricted to the interval q > 2k, while k/q < cos 9 < 1. This 
yields, after performing the ^-integral as shown in the appendix, 
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1 r°° Hf r°° 
A^h cp Jo v(ey J2k 
. f n k \f°(e nk )-f°(e nk + ahuj q )\ 
q 4 sinh (nu q /2kBT) 

where Io,o(q, 2/c/g) is given by Eq. (A. 5). We may carry out the final integration over energy 
by expanding the difference of the Fermi-functions in powers of uo q . Since the contribution 
due to terms involving higher order derivatives of the Fermi function is seen to vanish by 
the use of partial integration, we obtain 

1 f°° de dp 

([/ ' M) = «Vio 

(~ 2 q 2 + P)I , (q,2k/q)huj q 
2k q 4 sinh 2 (hw q /2k B T) ' 1 ' 

The remaining steps are standard. We assume that the temperature is much less than the 
Fermi temperature, so that the integral over e yields l/i>p, while k may be replaced by the 
Fermi momentum kp. Furthermore, we assume that the temperature is small compared to 
the characteristic temperature 9 defined by 

A; B = hk F c. (62) 

Then we obtain 

8 m 2 k F k B T ( 2 P \ 
^ pch 4 V" + 4fc|J 

Xe - 2fe |^/(i +7 ) e -20/r_ (63) 



(U,HU) 



This results in the final conductivity expression 



2e 2 

o = -^L, (64) 



where the length l a is given by 



, = jf'P^ e 2^/d + 7) e 2e/T - (65) 

m i (4fcp^, 2 + P) 

At low temperatures the conductivity thus increases exponentially, in agreement with the 
result of the numerical calculation (see Fig. 3). In the presence of impurities, as we shall 
see in the following subsection, the scattering against phonons yields a contribution to the 
inverse conductivity which is proportional to T 3 rather than exp(— 26/T), provided the 
temperature is sufficiently low that the impurities dominate the scattering. 
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C. Scattering from acoustical phonons and impurities 



When the impurity scattering dominates we may use the variational principle with a 
trial function which is proportional to the solution for impurity scattering alone. We shall 
consider the case where only the lowest Landau level, n = 0, is important. The trial function 
is thus chosen to be proportional to 1/ Bq+ as given by Eq. (|36|) , 

U{k) = k 2 exp (2k 2 l 2 h /[l + 7]) . (66) 

Since we have chosen the exact solution to the impurity problem as our trial function, the 
calculated upper bound on the contribution from electron-phonon scattering, l/cr p h, is exact, 
to lowest order in the magnitude of the electron-phonon coupling. At low temperatures we 
may make the approximation 

{U{k)-U{k')f^{^^ {k-k'f, (67) 

which is justified, since k — k! = q x and the restriction to low temperatures implies that 
the phonon momentum is small. Furthermore, we may neglect the deformation-potential 
coupling since the piezoelectric coupling dominates for small q according to Eq. fl53|), and 
also set the matrix element appearing in Eq. (|53|) equal to unity. By inserting U in Eq. ( |26|) 
and carrying out the integrals we obtain 

h ; _ 1 3C(3)l + 7 

74 




'^^Wl'^l 1^1 ■ (68) 



where Ep = ep/hu h , and 



Pm* 

£ = ( 69 ) 
pc z n 



This T 3 -behavior agrees well with our numerical calculations in the parameter range where 
it is expected to apply. In Figs. 3 and 4 we show the calculated conductivity, obtained by 
numerical solution of the integral equation, for samples with different amounts of impurities. 



VI. SCATTERING FROM OPTICAL PHONONS 

Next we investigate the combined effects of scattering by longitudinal optical phonons0 
and impurities. The chief differences from the previous section involve the g-dependence of 
the electron-phonon matrix-element and the absence of dispersion in the phonon frequencies. 
The conductivity is obtained by solving the Boltzmann equation in the limit where the 
temperature is much less than hu /k B . 
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A. Coupling to optical phonons 



We shall use the simple model where the phonon frequency is independent of momentum, 



u q ~ UJ . 

The electron-phonon matrix elemental is given by 

1 M 2 

\g(nn',q)\ 2 = - 2 ^M:'(u) 
where V is the normalization volume and where 

M 2 = 2vr^ 



(70) 



(71) 



(72) 



The function M is given by 



AC M 



[n 



k\e iq - r \n'k') 



(73) 



which is calculated in the appendix, where also u is defined. 

The approximation in Eq. (|70|) allows us to integrate over the y and z components of the 
phonon momentum. Using a delta-function for the wave-function in the ^-direction, we get 



J2 \9(nn',q) 



- K n (Qx), 



(2tt) 2 L 



(74) 



where 



dq y 



-°° \ Tx + QZ 



(75) 



B. Scattering from optical phonons and impurities 

We now consider the case of impurity and optical phonon scattering. The Boltzmann 
equation in this case reads 

e E , , , s (df nk \ mp fdf nk \° p 

The collision integral for the phonon scattering given in Eq. (|46|) simplifies in the case of 
optical phonons to 

(%0 OP = ~ w!ta ■ tf<* ,9k tS £ £ K'fa) - f(e nk + a^ ) 

V 9t /con ^27rL4smh (hw /2k B T) n , Qx a 

X (V'nfc - 4>n'k-q x ) ^n'fc-fe ~ ^nk ~ vTiUq), (77) 
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with K™'(q x ) given in Eq. (|7"5f) . The argument of the energy- conserving delta-function is 
zero for 

(78) 

for each value of nk, provided the square root is real. The branch index s assumes the values 
+1 and —1 corresponding to forward and backward scattering. The Boltzmann equation 
now simplifies to 

n' s' 

+ E <3?+ (e> ^(e) - s'Me + fftiuo)), (79) 



where the matrix B is given by Eq. (PS]). The optical-phonon scattering gives rise to the 
matrix C, which is given by 



■,1 r.' 



1 Mi 



" sV ' 1 /°(5)[l-r(£p 2 (27r)2 4sinh 2 (^o/2A:BT) 

x E \f°(e) - f{e + *H. )\ K\sk n (e) s>kAe + *H. )) _ 
a 1 v n {E)v n i{e + ahu Q ) 

In GaAs the typical optical phonon energies are about 36 meV. Furthermore, since the 
cyclotron energy %uo c in GaAs is 1.5 meV B/T, the typical situation will be that the energy 
spacing between the Landau levels as modified by the confinement potential, hut, is much 
smaller than the phonon energies, Huq 2> Two^. Since we are interested in temperatures that 
are small or comparable in magnitude to the energy spacing huh, we are therefore always 
in a situation where k^T is much less than Hujq. Under these circumstances it is possible to 
simplify the Boltzmann equation for electrical transport by considering it to be a coupled 
system of equations for the functions <p(e ± nhuo) with n integer. Since the current is 
mainly determined by the distribution of electrons <p(e) within a thermal layer of thickness 
comparable to ^b^, it is sufficiently accurate to neglect the contributions from <f>(e ±nhuj ), 
when n is different from zero. The collision integral may therefore be approximated by 
the "scattering-out" term, corresponding to the neglect of vertex-corrections in the Kubo- 
approach. 

In calculating the conductivity from Eq. ( [79] ) we thus neglect the "scattering-in" terms 
involving n /(e + ahooo), while retaining the full collision matrix for the impurities. We 
have then evaluated the solution of Eq. ( f79|) numerically and inserted the solution in the 
conductivity formula fllBp . Results are shown in Fig. 5. The low-temperature behavior is 
dominated by the impurities, which give rise to an initial increase of the conductivity, in 
accordance with Eq. (PH|). The optical phonons come in at higher temperatures, yielding a 
maximum of the conductivity as a function of temperature. The position of the maximum 
is roughly proportional to the logarithm of the strength of the impurity scattering, since 
the phonon contribution depends exponentially on temperature. Note however that the 
acoustical phonons have reduced the conductivity for the high mobility cases rather strongly 



16 



at the temperature range shown here, see Fig. 4. The dashed line represents the conductivity 
without optical phonon scattering. 

In comparing the curves Figs. 3-5 we observe that the optical phonon scattering in clean 
systems tends to dominate the total scattering already at fairly low temperatures, around 
50 K. This contrasts the situation in typical GaAs-based two-dimensional electron gases@, 
where the optical phonon scattering begins to dominate around 100 K. 



VII. CONCLUSION 



We have shown that the magnetoconductivity of quantum wires may be discussed in 
a simple and unified fashion within the framework of a Boltzmann equation, by taking 
into account the influence of the magnetic field on the electron group-velocity and the 
matrix-elements governing the transition probability. By treating in detail the scattering 
from acoustical as well as optical phonons in the presence of impurity scattering, we have 
determined the temperature dependence of the magnetoconductivity for realistic choices of 
parameters in GaAs-based structures. In particular, we have found that the magnetoconduc- 
tivity exhibits a maximum as a function of temperature, depending on the relative strength 
of the impurity and electron-phonon scattering. The calculated magnetoconductivity oscil- 
lates when the Fermi energy or the magnetic field is varied. Our detailed calculations show 
that the scattering against optical phonons in quantum wires is significant at temperatures 
somewhat smaller than the corresponding temperatures for the two-dimensional case. 

In order to relate our theoretical results to experimenting it is necessary to consider 
not only the diagonal conductivity element a xx , but also o xy . To compare with experiment 
one must therefore either determine a xy theoretically or compare our calculated a xx to 
simultaneous measurements of p xx and p xy . Alternatively, our results could be compared to 
experiments that use a Corbino-type geometry. The effects predicted in this paper should 
be observable in quantum wires of sufficient purity, since otherwise the electron density may 
vary considerably due to fluctuations in the electrostatic potential from the donors. 
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APPENDIX: 



In this appendix we study the square of the matrix element {nk\ exp(iq-r)\n'k'). Using 
the single electron wave functions in Eq. (Q) and performing the x- and ^-integrals one 
obtains 



\{nk\ exp(iq-r)\n k 



5k.k> 



n 



k,k'+q. 



mm- , u \n-ri\ 



X 



1 'max 

Ll n - n 'Hu) 



(Al) 
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where n m j n = min(n, n') and n max = max(n,n'); L™(u) are the Laguerre polynomials while 
u is given by 



^ = T —^Qj h r + MkY. (A2) 
1 + 7 I I 



The Kronecker delta appearing in Eq. ([Al]) , and which is written explicitly in the expression 
for the integral operator H in Eq. (f|4|), is a consequence of the translational invariance along 
the x-direction. 

In calculations involving acoustical phonons it is natural to use the polar coordinates 
(q,9,(p) with the x-axis as the polar axis so that q x = qcosO, q y = qsm9 coscj), and q z = 
q sin 6 sin 0. According to Eqs. ( f44[) and (|53[) the only ^-dependence is through u defined 
above. It is therefore of interest to calculate the integral 



I n , nl (q,co S 9) = £ ' d<Pu\ n - n '\ [L^\] 2 e~ u . (A3) 

From Eq. ( |A2|) it follows that the ^-dependence of u is of the form u = (3 + a cos 2 0, so the 
integrand of I n y therefore takes the form of a polynomial in cos 2 times exp(— a cos 2 4>). 
The integral over <fi may then be carried out explicitly in terms of the beta-function and 
Gauss' hypergeometric function. 

In the case where we consider only the lowest Landau level we just need to calculate J ,o- 
Using the results quoted above we find 

h,o{q,q x /q) = 2vrex P (^^(qxh) 2 -^-^ 2 h (q 2 -ql)j , 

xIoQW-a) (A4) 
where Iq(^) is the th Bessel function of imaginary argument. 
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TABLES 



TABLE I. The GaAs constants used in this paper. Unless otherwise indicated the values are 



taken from Ref. 27 



parameter symbol value 



ion mass density p 5.3x10 kg m 

longitudinal sound velocity c 5.2xl0 3 ms -1 

transverse sound velocity xc 3.0xl0 3 ms _1 

sound velocity ratio x 0.58 

static dielectric constant kq 12.8 

high-frequency dielectric constant Kqo 10.6 

piezoelectric constant! h\4 1.38xl0 9 Vm~ 1 

piezoelectric coupling-Eq. (|52"|) P 5.4x 10~ 20 J 2 m -2 

deformation potentialEil 3 2.2xlO~ 18 J 

effective electron mass m* 0.067 mo 

optical phonon energy %ujq 36 meV 
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FIGURES 



FIG. 1. The two lowest /c-dependent energy bands eofe an d £\k of the quantum wire. The dashed 
line represents an arbitrary Fermi level above the bottom of the second band. 



FIG. 2. The conductivity for a quantum wire with impurity scattering plotted versus the Fermi 
level for three different choices of the confinement potential corresponding to 7 = 1, 0.5, and 0.25. 
The dashed line is the T = result and the solid lines correspond to k^T = 0.05 hu>h- The magnetic 
field is 9 T. 



FIG. 3. Plots of the normalized conductivity versus temperature for a GaAs quantum wire at 
the magnetic field B = 9 T, the confinement parameter 7 = 0.5, and the Fermi level £f = 0.6 TioJh- 
Taking into account both impurity scattering, Eq. (j35|), and acoustical phonon scattering, Eqs. (46) 
and ([53]), the five dotted curves, a\, . . . , (T5, are numerical results for the zero- field, zero-temperature 
mobilities, /J.; m p = 7.5, 75, 300, 3000, and 30000 m 2 /(Vs) respectively. The last two rather unre- 
alistic high mobilities are considered to allow a study of the transition from impurity-dominated 
scattering to phonon-dominated scattering at temperatures low enough (T < ~ 4.0 K - see 
Eq. (|62l) ) for the approximative result Eq. (|6~4T ) to apply. The two full curves, <r a 4 and a a ^, are plots 
of the conductivity (rescaled to match 04 and 05) calculated from Eq. (|64|) where only the acousti- 
cal phonon scattering is present. The full curve o"i mp is the case where only impurity scattering is 
present. The two full curves, an and 0^5, are the results of assuming that the inverse conductivities 
for each scattering mechanism considered separately may be added, atj = cr a jai mp /(a a j + o"i mp ) 
(j = 4, 5), to approximate the exact numerical calculations of the total conductivity. 



FIG. 4. Plots of the normalized conductivity versus temperature for a GaAs quantum wire at 
the magnetic field B = 9 T and the Fermi level e-p = 0.6 huh- The confinement parameter in 
panel (a) is 7 = 1.0 and in (b) it is 7 = 0.5. The dashed curve in each panel is the case where 
only impurity scattering is taken into account. The four full curves in each panel is the result 
of combining impurity and acoustical phonon scattering for each of the following values of the 
zero-field, zero-temperature mobility, /ii mp = 0.9, 9, 90, and 900 m 2 /(Vs). 



FIG. 5. The conductivity with combined impurity and optical phonon scattering for a GaAs 
quantum wire. The normalized conductivity is plotted versus temperature for different choices of 
confinement potential strengths, Fermi energies and mobilities. The magnetic field is B = 9 T. 
The confinement parameter is 7 = 1 for panel (a) and (c) and 7 = 0.5 for (b) and (d); the Fermi 
level is ef = 0.6 huoh, i.e., close to the bottom of the first band, for (a) and (b), and panel (c) and 
(d) have £p = 1.3 fiuj^. The conductivity is shown for each of the following values of the zero-field, 
zero-temperature mobility, /^imp = 

0.9, 9, 90, and 900 m 2 /(Vs). The dashed curve in each panel is 
the case where only impurity scattering is taken into account. 



21 



